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Resumo 



Foram feitas simulagSes baseadas no modelo Penna para envelhecimento biologico, mas com o 
objetivo de estudar a especiagao simpatrica, isto e, a divisao de uma especie em duas ou mais 
populag5es, isoladas reprodutivamente, mas sem que nenhuma barreira ffsica as separe. Para 
tal, introduzimos urn novo tipo de competigao entre os individuos, atraves de um fator de Ver- 
hulst modificado. Esta nova competigao depende de certas caracteristicas fenotfpicas de cada 
individuo, que sao representadas por pares de tiras de bits. Estas tiras de bits sao lidas em par- 
alelo e nao tem estrutura de idade. Desta forma, o genoma de cada indivfduo tem duas partes: 
uma estruturada por idade e relacionada ao aparecimento de doengas hereditarias e outra, sem 
estrutura, ligada a competigao por alimento. Tambem foi introduzida uma selegao sexual no 
modelo, representada por um par de tiras de bits nao estruturadas independente. 

Nesta tese apresentamos tres modelos diferentes; dois deles se utilizam, alem da competigao 
entre os individuos, de uma ecologia que muda abruptamente, para obter a especiagao. Foram 
inspirados na especiagao dos tentilhSes de Darwin, uma familia de passaros que habitant as 
Ilhas Galapagos e nos ciclfdeos, uma familia de peixes que vivem nos Lagos da Nicaragua e 
no Lago Vitoria, na Africa. O terceiro modelo nao utiliza nenhuma mudanga ecologica: a 
especiagao e obtida dependendo apenas do grau de competigao entre grupos de individuos com 
caracteristicas fenotfpicas semelhantes. 

No modelo para descrever o processo de especiagao dos tentilhoes utilizamos apenas um 
par de tiras de bits para descrever uma unica caracteristica fenotfpica, que esta ligada tanto 
a selegao ecologica como a selegao sexual. Esta caracteristica equivale a morfologia dos bi- 
cos. Estudamos entao a mudanga na distribuigao destas morfologias em fungao das regras de 
acasalamento, dos valores dos parametros do modelo e da posigao da especie na cadeia ali- 



2 



mentar. Neste caso encontramos comportamentos qualitativos equivalentes aos que vein sendo 
estudados nestes passaros. 

Para simular os cicKdeos, utilizamos dois pares de tiras de bits independentes; um deles 
representa o tipo de maxilar do peixe e esta diretamente ligado a competigao por alimento. O 
outro par representa a cor, que e o trago ligado a selegao sexual. O objetivo principal foi o 
de analisar a estabilidade das distribuigSes tanto dos maxilares como das cores destes peixes. 
Obtivemos a especiagao simpatrica com este modelo, embora nao exatamente aquela observada 
nos ciclideos. Verificamos ainda que se estabelece uma forte correlagao entre os dois tragos, 
sempre que a especiagao ocorre. 

Com o modelo que estuda a especiagao simpatrica em fungao do grau de competigao numa 
ecologia fixa, encontramos uma transigao de fase que possibilitara o estudo da reversibilidade 
ou nao deste processo de especiagao. 
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Abstract 



We perform simulations based on the Penna model for biological ageing, now with the purpose 
of studying sympatric speciation, that is, the division of a single species into two or more 
populations, reproductively isolated, but without any physical barrier separating them. For that 
we introduce a new kind of competition among the individuals, using a modified Verhulst factor. 
The new competition depends on some specific phenotypic characteristic of each individual, 
which is represented by a pair of bitstrings. These strings are read in parallel and have no 
age structure. In this way, each individual genome consists of two parts. The first one has 
an age-structure and is related to the appearance of inherited diseases; the second part is not 
structured and takes into account the competition for the available resources. We also introduce 
sexual selection into the model, making use of another non-structured and independent pair of 
bitstrings. 

In this thesis we present three different models; two of them use, besides the competition, 
a sudden change in the ecology to obtain speciation. They were motivated by the speciation 
process observed in the Darwin finches, a family of birds that inhabits the Galapagos Islands, 
and also by that observed in the cichlids, a family of fish that lives in the Nicaragua Lakes and 
in the Vitoria Lake, in Africa. The third model does not use any ecological change: sympatric 
speciation is obtained depending only on the strength of competition among individuals with 
similar phenotypic characteristics. 

In the model to describe the speciation process of the finches we use a single pair of non 
age-structured bitstrings to describe one phenotypic characteristic which is related to both 
competition and mating preference. This characteristic corresponds to the beak morphology. 
We study the changes in the distribution of the beaks according to the mating rules, to the 
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values of the parameters and to the position of the species in a food-chain. In this case we obtain 
qualitative results that are in agreement with field observations of these birds. To simulate the 
cichlids, we use two pairs of non-structured and independent bitstrings. One pair represents 
the fish type of jaw and is directly related to competition for food. The other pair represents 
the colour, which is the trait related to sexual selection. Our main purpose has been to analyze 
the stability of the distributions of both jaws and colours. In this case we obtain sympatric 
spectiation, although not exactly the one experimentally observed in the cichlids. 

Using the model without ecological changes, we study the probability of speciation as a 
function of the competition strength and find a phase transition, which may allow the study of 
the reversibility or not of the speciation process. 
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Introdugao 



processo de especiagao simpatrica, cuja definigao sera vista mais adiante, esta baseado na 
teoria da evolugao de Darwin. Esta teoria tem uma matematica pouco formalizada, nao so 
por tratar de sistemas adaptativos mas tambem devido ao tipo de interagao local entre os 
individuos, que inclui o cruzamento e a recombinagao de seus genomas. Erwin Schroedinger, 
no seu livro What is life? 0, foi urn dos primeiros cientistas a tentar interpretar e conjecturar 
a vida em termos dos conhecimentos da quhnica e da fisica. 




Figura 1: Erwin Schroedinger escreveu, em 1948, o primeiro livro dedicado a problemas de 
Biologia sob o ponto de vista fisico. 

Numa palestra a respeito de diferentes aspectos das celulas vivas, Schroedinger explica a 
uma plateia de leigos que o assunto era dificil, embora a dedugao matematica pouco fosse 
utilizada: "A razao disto nao e que o assunto seja simples o bastante para poder ser explica- 
do sem matematica, mas sim porque e complexo demais para ser completamente accessi'vel a 
matematica." Atualmente um grande numero de fisicos tem investido no assunto, principal- 
mente utilizando computadores, os quais funcionam como um grande laboratorio no qual os 
passos da evolugao podem ser recriados e medidas estatisticas podem ser feitas. 
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Evolucao das especies 

O interesse pela evolugao organica teve suas origens entre os Gregos e Romanos. As ideias 
de Empedocles (cerca de 490-435 a.C.) [2], por exemplo, ja demonstravam serios intentos de 
compreender as mudangas biologicas. No entanto, ate meados do seculo XIX muitos obstaculos 
se interpuseram no caminho de uma apreciagao cientffica racional da evolugao. Um dos mais 
importantes obstaculos foi a crenga, amplamente aceita, da "imutabilidade das especies", que 
considerava todo tipo de formas de vida como entidades estaticas, sem mudangas no passado e 
sem possibilidades de mudanga no futuro. A obtengao sistematica de provas e a construgao de 
uma teoria de evolugao das especies por selegao natural foi conseguida por Charles Darwin, como 
resultado de 20 anos de trabalho arduo. Ela foi apresentada em 1857, ao mesmo tempo em que 
Alfred Russel Wallace tambem o fazia, ja que chegara as mesmas conclusoes de Darwin. Ambos 
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Figura 2: Representagdo esquemdtica dos principals argumentos em favor da evolugao por 
selegao natural defendidos por Charles Darwin e Alfred Russell Wallace; tabela de Wallace. 

argumentaram que mesmo que existissem especies com morfologias fixas, estas haviam aparecido 
por selegao natural entre os membros variaveis de especies anteriores. Um interessante esquema 
dos conceitos por eles desenvolvidos pode ser visto na figuraEl Nessa epoca nao se sabia explicar 
como as variagoes podiam ser transmitidas aos descendentes, pois os estudos desenvolvidos 
por Gregorio Mendel (1866) nao tinham tido repercusao. Somente em 1900, 34 anos depois, 
tres botanicos redescobriram as mesmas leis basicas da her edit ariedade apontadas por Mendel 
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jHj. A adapatagao e a selegao natural, embora de facil comprensao intuitiva, sao na verdade 
bastante complexas. Estuda-las significa compreender nao apenas relacionamentos ecologicos 
intrincados, como tambem a matematica avangada envolvida na genetica das populagSes |UEj- 
Atualmente entendemos que a evolugao tern como base a mudanga na constituigao genetica 
da populagao. E sobre estas diferengas geneticas, originadas por mutagSes e recombinagao, que 
atua a selegao natural j2]. Os genes, tal como sao entendidos hoje em dia, sao segoes de 
cadeias de DNA responsaveis pelo funcionamento e aparencia dos organismos. A distingao 
entre os efeitos dos genes num organismo e os genes propriamente ditos e considerada atraves 
dos termos fenotipo e genotipo, respectivamente (Futuyma 1979) jE]. O genotipo content as 
informag5es geneticas, que ao serem "compiladas" , dao origem aos fenotipos. 

Especiacao 




Figura 3: O Istmo do Panama isolou uma especie do genero de lagostas Alpheus, ha tres milhoes 
de anos atrds, sendo um dos exemplos mais documentado de especiacao alopdtrica. 

Especiagao e um processo pelo qual uma populagao original se divide em duas ou mais, sem 
que haja um fluxo de genes entre elas (sexualmente isoladas). Os dois tipos mais importantes 
de especiagao sao a alopatrica e a simpatrica. A especiagao alopatrica basicamente consiste na 
separagao geograflca da populagao em duas ou mais partes. Esta separagao pode ser devida 
a migragao de uma das partes ou ao aparecimento de um obstaculo ffsico que impega a troca 
de gens entre elas. A diferenciagao destas novas especies e bem conhecida, e deve-se a fatores 
estocasticos (deriva genetica) e a processos seletivos de adapatagao. O melhor exemplo docu- 
mentado deste tipo de especiagao foi aquele produzido na formagao do Istmo do Panama, ha 
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tres milhoes de anos atras. O aparecimento desta barreira geografica separou as populagoes de 
organismos aquaticos entre as aguas dos oceanos Pacifico e Atlantico. Dentro do genero Alpheus 
de lagostas, por exemplo, deu origem a 7 especies distintas, com representantes de cada lado 
do Istmo. Estas especies sao consideradas gemeas, pois elas diferem apenas ligeramente na 
morfologia, mas estao isoladas geneticamente. 

O segundo tipo de especiagao e a chamada especiagao simpatrica, que corresponde a divisao 
da populagao dentro de urn mesmo espago geografico. Esta divisao e produzida pela selegao 
natural, que atua atraves da competigao pelas fontes de alimento em conjunto com a selegao 
sexual. A busca pelo entendimento deste tipo de especiagao tern gerado muitos modelos teoricos 
e computacionais; no nosso caso, estes ultimos correspondem a ferramenta que temos utilizado 
para estudar este fascinante processo. 

Os modelos computacionais usados nos capftulos dois, tres e quatro, tern uma origem co- 
mum, que esta descrita no primeiro capitulo. Com o tipo de competigao que utilizamos obser- 
vamos que existe um grau de competigao critico a partir do qual ocorre ou nao a especiagao. Co- 
mo sera visto no capitulo 4, para graus menores que o critico nao e possivel achar especiagao 
simpatrica e para valores maiores, sempre e possivel encontra-la. Nos capftulos 2 e 3 trabal- 
hamos exatamente com este grau critico de competigao, onde a especiagao pode ou nao ocorrer, 
devido as flutuagSes nas densidades populacionais. Afinal, e possivel encontrar especiagao 
simpatrica na natureza, mas nao o tempo todo. As simulag5es feitas foram inspiradas nos ex- 
emplos dados a seguir, observados experimentalmente. Maiores detalhes sobre estes processos 
de especiagao serao vistos nos capftulos correspondentes. 

Exemplos de especiagao 

Um dos mais famosos exemplos de especiagao sao os classicos tentilhSes (finches) de Darwin, 
das ilhas Galapagos, que constituem um grupo de 14 especies com formas e tamanhos de 
bicos diferentes. O casal Peter Grant e Rosemary Grant, da Universidade de Princeton, EUA, 
comegou em 1973 a estudar a populagao de duas especies de uma pequena ilha, a Daphne 
Maior IB] • Eles contaram as populagoes e mediram certos tragos do tentilhao terrestre de 
bico medio (Geospiza fortis) e do tentilhao dos cactos (Geospiza scandens). Os tragos basicos 
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eram o tamanho do corpo e o tamanho e a forma do bico. Eles notaram uma forte correlagao 
dessas medidas com eventos naturais pelos quais a ilha passou. Ao fazer o estudo de longo 
prazo, foi possfvel medir o impacto nos tentilhSes dos momentos de seca intensa ou de chuvas 
torrenciais - como a provocada pelo fenomeno climatico El Nino, em 1983. No caso do G. fortis, 
por exemplo, notaram que o tamanho medio do bico aumentava nos anos de seca, ver figura 
quando apenas sementes maiores e duras estavam dispom'veis. Em tempos mais limidos, bicos 
menores eram mais comuns. Eles tambem descobriram que estes passaros acasalam atraves do 
tipo de canto que os machos emitem para atrair as femeas 
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Figura 4: Notdvel mudanga no tamanho dos bicos da especie Geospiza em fungao das mudangas 
climdticas. 

Como todo bom exemplo classico, os tentilh5es tern ajudado a entender os efeitos da se- 
legao ecologica, devida a competigao por comida, e da selegao sexual, durante um processo de 
especiagao. As simulagoes inspiradas nestes passaros podem ser vistas no capftulo 2. 

O capftulo 3 corresponde a especiagao da famflia de peixes chamados ciclfdeos, figura E] 
Existe uma forte evidencia, que o tipo de processo de especiagao destes peixes e simpatrico 
[TUj . Contudo, como sera visto no capftulo 3, as origens deste processo nos cichdeos tern 
dividido em dois grupos aqueles que estudam especiagao simpatrica. O primeiro grupo acha 
que a especiagao simpatrica e guiada pela divergencia causada pela selegao natural devida a 
competigao pelas fontes de alimento. O outro grupo defende que este tipo de especiagao e 
guiado pela divergencia causada por selegao sexual. 
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Figura 5: Dados encontrados no trabalho de Wilson [11]. Lado esquerdo: Mapa da Nicaragua 
que mostra o lugar dos lagos onde foram colhidas as amostras da jamilia de peixes chamados 
ciclideos. Lado direito: A variagao da morfologia de cores e maxilares destes peixes, encontrados 
nos lagos da Nicaragua. 




Figura 6: A larva da raosca da maga e um outro exemplo recente de processo de especiagao 
simpdtrica em curso. Lado esquerdo: A mosca da fruta Rhagoletis pomonella adulta. Lado 
direito: A larva da mosca da maga; dependendo da raga da mosca, a larva pode ser deixada na 
superficie ou no centro da maga. 
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Finalmente o capftulo 5 refere-se ao estudo da obtengao da especiagao simpatrica numa 
ecologia fixa, variando-se apenas o grau de competigao entre os indivfduos. Encontramos uma 
transigao de fase cuja ordem, ainda nao identificada, pode ser bastante relevante no que diz 
respeito a reversibilidade ou nao do processo de especiagao. Neste capftulo a motivagao bi- 
ologica veio da mosca da maga - ver figura |3 O processo de especiagao simpatrica em curso 
destas moscas parece ser reversfvel, mas a bibliografia a respeito do mesmo ainda esta sendo 
acumulada, pois sua observagao e bastante recente. 

A tese e conclmda com observagoes gerais dos resultados das simulagoes apresentadas nos 
capitulos 2,3 e 4. 
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Capitulo 1 
Modelos 



Conforme mencionado na introdugao, os ingredientes principals para se obter a especiagao 
simpatrica sao a disputa por alimento e a selegao sexual. Ambos os ingredientes estao rela- 
cionados a caracteristicas fenotipicas externas do individuo, tais como tamanho e cor. Nossas 
simulagSes sobre especiagao simpatrica baseam-se na introdugao de uma nova parte ao "genoma 
cronologico" que representa o individuo no modelo Penna , a fi m de atribuir ao mesmo um 
"fenotipo". 

O modelo Penna e baseado na teoria da selegao natural de Darwin para a evolugao das 
especies e na teoria do acumulo de mutagoes para explicar o envelhecimento biologico. Desde 
sua publicagao em 1995, vem sendo utilizado com sucesso na compreensao de muitos fenomenos 
evolucionarios observados na natureza, tais como a senescencia catastrofica do salmao, a auto- 
organizagao da menopausa, as vantagens da reprodugao sexuada, etc... Para uma revisao das 
diversas aplicagoes teoricas do modelo Penna veja [T3J E] ; para possiveis aplicag5es praticas 
veja [13 e para solugoes analiticas veja [T6| IT7j. 

Embora o estado final do modelo Penna, calculado analiticamente ^Bj, seja a morte da 
populagao toda, o tempo que este estado final leva para ser atingido cresce exponencialmente 
com o tamanho da populagao. Assim, os individuos vivem num estado quase-estacionario de 
duragao quase infinita, onde a populagao apresenta uma distribuigao de individuos por idade 
que se mantem constante durante todas as simulagSes, por mais longas que estas sejam (bilh5es 
de passos). Esta propriedade do modelo certamente facilita o estudo da dinamica de populag5es 
em processo de especiagao. 
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1.1 O modelo Penna 



Na versao sexuada do modelo Penna o genoma de cada individuo e representado por duas tiras 
de 32 bits cada, que sao lidas em paralelo. Elas content a informagao de quando os sintomas de 
uma dada doenga hereditaria vao aparecer, sendo por isto chamadas de "genoma cronologico" . 
Cada uma das tiras content a heranga genetica de um dos pais, sendo as doengas representadas 
por bits 1. Se um dado individuo possui dois bits iguais a 1, por exemplo na terceira posigao de 
ambas as tiras (homozigoto), isto indica que aquele individuo vai comegar a sofrer dos sintomas 
de uma doenga no terceiro perfodo de sua vida. Assim, cada individuo pode viver no maximo 
por 32 periodos. Se o individuo for heterozigoto numa dada posigao (bit 1 numa tira e zero na 
outra), ele so ficara doente se naquela posigao o bit 1 for dominante. No inicio da simulagao 
defme-se quantas posigoes serao dominantes e sorteia-se aleatoriamente quais serao elas. Estas 
posig5es sao as mesmas para todos os genomas e sao mantidas fixas durante todo o processo 
de evolugao da populagao. Um passo computacional significa ler mais um bit do genoma de 
todos os individuos. Se em qualquer passo, o mimero de doengas acumuladas num dado genoma 
atinge o limite L, aquele individuo morre. 

A femea que sobrevive ate a idade minima de reprodugao R acasalar aleatoriamente 

e a cada passo, com um macho de idade tambem > R, gerando NF filhos a cada vez. O genoma 
do filho e construido atraves do cruzamento e da recombinagao das tiras dos pais, conforme 
esquematizado na figura 11.11 Primeiro corta-se o genoma da mae numa posigao aleatoria e 
une-se dois pedagos complementares para formar o gameta feminino (uma unica tira de bits). 
A seguir, introduz-se aleatoriamente M mutag5es nocivas, 1 cujas posigoes estao indicadas por 
setas na figura [TTT1 Neste caso, se o bit escolhido para mutar no gameta feminino vale 0, ele se 
torna 1 no genoma do filho, mas se o bit escolhido ja for 1, ele permanece inalterado no filho 
(nao ocorre mutagao). O mesmo processo de cruzamento e recombinagao com taxa de mutagao 
M ocorre no genoma do pai, gerando o gameta masculino. O genoma do filho e o resultado da 

uniao dos dois gametas e o sexo e aleatoriamente escolhido, com 50% de chance para cada um. 

1 As mutagoes normalmente ocorrem durante a duplicagao do DNA e devem-se a perda ou adigao de uma 
base, ou a substituigao de uma base por outra. A maioria das mutagoes sao neutras, mas as mutagoes ruins 
sao 100 vezes mais provaveis que as boas (corretoras) |18| . sendo que estas ultimas em geral dao origem a uma 
nova especie. 
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Figura 1.1: Formagdo dos gametas feminino e masculino que ddo origem ao genoma do filho. 
O processo de cruzamento e recombinagdo (lado esquerdo ) e feito da mesma forma tanto no 
genoma do pai quanto no da mae. As setas indicam onde mutagoes nocivas aleatorias ocor- 
reram. O desenho (lado direito) mostra que uns dos elementos principals para a continuidade 
da vida e o processo de divisao celular chamado de meiose, que ocorre nas celulas germinativas 
e e o unico onde acontece um cruzamento entre dois cromossomos homologos. 



Mesmo permitindo apenas mutagoes ruins, a populagao que se obtem com a dinamica ate 
aqui descrita cresce exponencialmente. Para evitar tal comportamento explosivo, a versao origi- 
nal do modelo Penna adota o chamado fator logistico de Verhulst, dado por V(t) = P(t) /MAX, 
onde P(t) e a populagao total no tempo t e MAX e a capacidade de sustentagao do ambiente. 
Assim, em cada passo e para cada individuo gera-se um numero aleatorio entre e 1, cujo valor 
e comparado com o valor de V. Se o numero gerado for menor que V, o individuo morre, inde- 
pendente da idade ou do genoma. Este termo tambem representa uma competigao aleatoria por 
espago e comida, que depende do numero maximo, MAX, de individuos que o ambiente pode 
sustentar. Exatamente por ser um termo totalmente aleatorio, o fator de Verhulst vem sendo 
criticado na literatura [19 . A introdugao de um fenotipo no modelo Penna evita a utilizagao 
de tal termo, como veremos a seguir. 
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1.2 Modelos com fenotipo. 



Para simular uma caracteristica fenotipica do individuo adotamos o mesmo procedimento intro- 
duzido por Medeiros, Sa Martins, S. Moss de Oliveira e Medeiros j20], isto e, acrescentamos mais 
urn par de tiras de bits ao genoma cronologico do modelo Penna. Este novo par nao apresenta 
estrutura de idade, mas suas tiras tambem sao lidas em paralelo. No momento da reprodugao 
este par sofre o mesmo processo de cruzamento e recombinagao indicado na figura lITTl e existem 
tambem posig5es dominantes e recessivas para os bits 1. Este numero de posigoes dominantes e 
dominado D, para parte com estrutura e DF para parte sem estrutura. A unica diferenga esta 
nas mutagoes, pois para este novo par elas podem se dar nas duas diregoes (0 <-> 1), com uma 
taxa MF por bit (dos 32 bits escolhe-se so 1 bit). O numero total k (g [0, 32]) de bits iguais a 1 
e contado nestas duas tiras, levando-se em conta as posigoes dominantes. Este numero e entao 
relacionado a algum tipo de caraterfstica fenotipica do individuo, conhecida na biologia como 
trago. Supondo que tal caracteristica seja, por exemplo, o tamanho do animal ou o tamanho da 
cauda de um passaro, mutagSes para um lado ou para o outro apenas aumentam ou diminuem 
tal tamanho nas proximas gerag5es, o que em principio nao pode ser considerado bom nem 
ruim. O carater seletivo destas mutagoes so sera defmido quando for estabelecida a ecologia a 
qual a populagao estara sujeita e/ou a preferencia sexual por determinado fenotipo na hora do 
acasalamento. A figura fOl ilustra as duas partes do genoma de cada individuo e como e feita 
a leitura de cada parte. 

Ao nos referirmos apenas a parte nao estruturada do modelo como fenotipo, estamos come- 
tendo um abuso de linguagem. Como mencionado na introdugao, qualquer tradugao ou compi- 
lagao do genotipo corresponde a um fenotipo. Assim, a leitura da parte estruturada por idade 
do modelo Penna e representada por uma tira de bits na figura ll.2[ lado esquerdo, tambem 
corresponde a um fenotipo. Contudo, para distinguir com facilidade a parte estruturada da nao 
estruturada, chamaremos de fenotipo apenas a tradugao da parte nao estruturada do genoma. 

Na natureza existe uma enorme variedade de mecanismos relacionados a capacidade do 
individuo com um dado fenotipo de produzir uma descendencia apta a sobreviver. Um destes 
mecanismos e o da selegao natural na competigao por diferentes fontes de alimento, que atua no 
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Figura 1.2: Representagao das partes cronoldgica e sem estrutura de idade dos genomas. As 
setas indicam as posigoes onde os bits 1 sao dominantes e a leitura do genoma e feita de acordo 
com tais posigoes. Para a parte estruturada por idade, tanto o numero de bits 1 quanto suas 
respectivas posigoes sao relevantes; jd para a parte nao estruturada, apenas o numero efetivo de 
bits 1 importa. O desenho do lado direito mostra uma caricatura do pai da Genetica Moderna, 
Gregor Mendel, que introduziu os conceitos hoje conhecidos como domindncia e recessividade 
nos estudos estaMsticos que realizou com ervilhas em 1865. 

indivfduo de acordo com seu trago ecologico, isto e, com sua caracteristica fenotipica relevante 
para a obtengao de alimentos. Outro mecanismo e o da selegao sexual, no indiviiduo 

de acordo com alguma caracteristica fenotipica relativa a preferencia sexual no momento do 
acasalamento - o trago sexual. A selegao sexual e um ingrediente fundamental para se obter 
o isolamento reprodutivo na especiagao simpatrica. Para as especies cujo trago ecologico e o 
mesmo que o sexual, nos faremos uso de apenas um par de tiras de bits nao estruturadas por 
idade. No caso em que a selegao natural devido a ecologia atua num trago que e diferente 
daquele utilizado para definir a preferencia sexual, serao usados dois pares independentes de 
tiras de bits nao estruturadas. 

1.2.1 O trago ecologico 

A caracteristica fenotipica de cada indivfduo ligada a competigao pelo alimento disponivel fica 
entao defmida pelo numero k de bits 1 efetivos no segundo par de tiras de bits de seu genoma 
(0 < k < 32). Em nossas simulagoes dividimos a populagao em tres grupos, de acordo com 
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os valores de k. O grupo Pi consiste dos individuos com < k < n\\ o grupo do meio, P m , 
corresponde aos fenotipos intermediarios onde rii < k < n 2 e o terceiro grupo, P 2 , corresponde 
aos fenotipos com n 2 < k < 32. Se supusermos que o valor de k esta ligado ao tamanho dos 
indivfduos, os grupos Pi e P 2 correspondem aos tamanhos extremos, sendo os individuos do 
grupo Pi bem pequenos (com poucos bits 1 efetivos) e os do grupo P 2 bem grandes (ricos em 
bits 1). 

O carater seletivo deste trago e dado pelo seguinte fator de Verhulst modificado: 

y Niimero de individuos com fenotipos que competem pela mesma fonte de comida 

grupo ' Capacidade de sustentagao da fonte de alimento dispomvel F(k,t) 

A seguinte competigao intra-especifica e mantida fixa em todas as simulagoes, mas a ca- 
pacidade de sustentagao da fonte de alimento dispomvel, F(k, t), pode ou nao variar no tempo. 
Com este novo fator de Verhulst as mortes por disputa de comida deixam de ser totalmente 
aleatorias e passam a depender dos fenotipos e da fonte de alimento dispomvel. 

A forma explicita da competigao inter-especifica, entre os individuos de cada grupo e dada 
por 12] : 



Vi(k,t) se < k < Hi, 

V(k,t) = { V m {k,t) seni<k<n 2l (1.1) 

V 2 (k,t) sen 2 <k<32. 

com a competigao intra-especifica igual a: 



Vi(k,t) 



Pi(k,t) + P m (k,t) 
F(k,t) 



v (h f s P m {k,t)+X * jPi(k,t) + P 2 (k,t)] 

Vm^t) - , [1.2) 



V 2 (k,t) 



p 2 (k,t) + p m (k,t) 



F(k,t) 

onde Pi(k, t) e a populagao com fenotipos em um dos extremos, P m (k, t) e a populagao do meio 
e P 2 (k,t) e a populagao com os fenotipos no outro extreme O parametro X regula com que 
fragao das populagoes dos extremos os individuos de fenotipos intermediarios irao competir. Ja 
os individuos do grupo Pi competem entre si e com todos aqueles de fenotipos intermediarios, 
mas nao competem com os fenotipos do outro extremo; o mesmo ocorre com os individuos 
do grupo P2, que tambem nao competem com os do grupo Pi. Em cada passo e para cada 
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indivfduo com fenotipo k gera-se urn numero aleatorio entre e 1, cujo valor e comparado com 
o valor de V(k,t). Se o numero gerado for menor que V(k,t), o indivfduo morre. 

Este tipo de competigao ja e bastante adotado na literatura biologica da especiagao simpatrica 
[2] , porque produz um aumento na diferenciagao genetica entre grupos com fenotipos extremos, 
abrindo o caminho para o isolamento reprodutivo nestes grupos quando a selegao sexual e 
introduzida. 

1.2.2 O trago sexual 

A selegao sexual e definida atraves do acasalamento nao aleatorio. Adicionamos ao genoma de 
cada femea mais um gene, representado por um unico bit, que determina se a femea e seletiva 
ou nao no momento da reprodugao. Se este bit for igual a 1, ela escolhe o parceiro de acordo 
com alguma preferencia relacionada ao trago sexual; caso contrario, ela escolhe aleatoriamente 
um macho para acasalar. Este gene e herdado da mae pela filha com uma probabilidade de 
mutagao (em ambas as diregoes) igual a MS. Isto significa que as filhas de femeas seletivas 
podem ou nao ser seletivas e vice-versa. As tres regras de selegao sexual que utilizaremos para 
obter o isolamento reprodutivo sao: 

• Regra com forte diregao : As femeas seletivas com fenotipos extremos escolhem entre 
um numero NM de machos, aquele com o fenotipo mais extremo do mesmo grupo que 
o dela. Supondo que o trago sexual seja tambem o tamanho do animal, por exemplo, 
uma femea grande tenta acasalar com o maior dos machos do conjunto NM, ao passo 
que uma femea pequena escolhe o menor dos machos do conjunto. As femeas seletivas 
com fenotipos intermediaries atuam como as femeas seletivas de um dos dois conjuntos 
de fenotipos extremos, com 50% de probabilidade para cada conjunto. 

• Regra com diregao: Femeas seletivas com fenotipos extremos acasalam com o primeiro 
macho sorteado aleatoriamente que pertenga ao mesmo grupo fenotipico que ela. As 
femeas seletivas com fenotipos intermediaries se comportam como aquelas seletivas de 
um dos conjuntos extremos, defmido a priori. 
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• Regra da diferenga : A femea seletiva com fenotipo kp escolhe, dentre um numero NM 
de machos, aquele cujo fenotipo ku apresenta a menor diferenga para o seu proprio kp, 
isto e, 

mmimo{|/c F - k Ml \, \k F - k Ma \, ■•■,\k F - k MNM \} 

Se o trago sexual for diferente do ecologico, acrescentamos mais um par de tiras de bits nao 
estruturadas por idade ao genoma do individuo, cujo processo de cruzamento, recombinagao e 
mutagao e o mesmo do trago ecologico. 

A nomenclatura utilizada ate aqui para descrever os modelos sera mantida nos capitulos 
seguintes. 

Em todas as simulagoes a populagao inicial, de 3000 femeas e 3000 machos, tem a parte 
cronologica do genoma constituida apenas de zeros, isto e, todos os individuos sao bebes livres 
de doengas hereditarias. A parte sem estrutura de idade, que define as caracteristicas externas 
dos individuos ligadas a ecologia e a preferencia sexual e constituida de zeros e uns, aleato- 
riamente distribuidos. Para garantir que nossas medidas fossem tomadas em populag5es ja 
equilibradas geneticamente, acompanhavamos sempre a distribuigao de idades dos individuos e 
so computavamos os dados muitos passos depois da mesma ter se tornado constante. 
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Capitulo 2 

Modelos com um trago fenotipico 



Com o tipo de dinanima descrita no capitulo anterior e inspirados num processo de especiagao 
observado nos tentilhoes das ilhas Galapagos, nos estudamos os possfveis mecanismos de espe- 
ciagao simpatrica em cadeias alimentares, utilizando populagoes com um linico trago fenotipico 
e vivendo numa ecologia onde a fonte de alimento muda com o tempo. 

Primeiro analisamos detalhadamente o processo de especiagao simpatrica de uma populagao 
isolada de herbivoros, que sao diretamente afetados pelas mudangas morfologicas que ocorrem 
na fonte de alimento, causadas, por exemplo, por oscilagSes sazonais no regime de chuvas. 
Tentando ficar o mais perto possfvel da realidade e levando em conta que na natureza existem 
predadores que comem herbivoros, os quais por sua vez comem plantas, nos comparamos os 
processo de especiagao simpatrica em cadeias alimentares de duas e de tres especies. 

2.1 Motivacao biologica 

Existem atualmente varios trabalhos experimentais e teoricos sendo desenvolvidos acerca da 
evolugao dos tentilhoes, tambem conhecidos como passaros de Darwin. A significativa bibli- 
ografia das observagSes feitas nas ilhas Galapagos conta com estatisticas do tamanho do corpo, 
morfologia do bico e series temporais dos sons destes passaros (ver figuraEJ). As medidas 
mostram que a diversificagao na morfologia dos bicos tern formado padroes de evolugao vocal 
que indicam uma correlagao entre o tamanho do bico e o tipo de canto [22 123] • Nesses trabal- 
hos tambem e mencionado que a restrigao fisica que a morfologia do bico causa na emissao do 
tipo de som dos tentilhoes pode ter levado a populagao original a um isolamento reprodutivo e 
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consequentemente, a uma rapida especiagao [21] . 

Nas simulagoes a seguir modelamos a especiagao simpatrica em populagoes onde existe uma 
relagao direta entre "o tamanho do bico" e "o som emitido" pelos "passaros" . Isto e, populag5es 
cujo trago ecologico e o mesmo que o sexual (tamanho do bico) e portanto necessitamos apenas 
de um par de tiras de bits nao estruturadas por idade. 

2.2 Especiagao para herbfvoros 

Inicialmente estudaremos uma populagao de passaros herbfvoros que se alimentam de plantas ou 
sementes cuja morfologia varia no tempo. Primeiro analisaremos a competigao intra-especifica 
pelas fontes de alimento. A seguir introduziremos a selegao sexual que, em conjunto com a 
disputa por alimento, levara a populagao a especiar. Isto e, a populagao inicial de fenotipo 
intermediario se dividira em duas populagoes com fenotipos extremos diferentes e sexualmente 
isoladas. 

2.2.1 Competigao sem selegao sexual 

Utilizaremos uma ecologia na qual a fonte de alimento muda de uma fonte que e abundante 
para todos os fenotipos para outra na qual os fenotipos intermediaries sao desfavorecidos em 
relagao aos extremos. A capacidade de sustentagao desta fonte varia no tempo da seguinte 
forma: 

[ 1 V < t < 4 x 10 3 , 

F(k,t) = 2 x 10 5 x J (2.1) 
[ 0.1+^|=^ V4x 10 3 <t<3x 10 4 , 

onde /ceo fenotipo da populagao de herbfvoros e t e o passo da simulagao. 

Os fatores de Verhulst aos quais os individuos estao sujeitos sao dados pela equacao lTT2l Para 
t < 4 x 10 3 a capacidade de sustentagao do ambiente independe dos fenotipos. Corresponde, por 
exemplo, a uma epoca de chuvas onde existe uma grande variedade de tamanhos de plantas ou 
sementes. A partir do passo 4 x 10 3 a capacidade de sustentagao passa a depender dos fenotipos, 
ou melhor, cada passaro passa a ter maior ou menor facilidade de arrumar comida, dependendo 
de seu valor de k. Este periodo corresponderia as epocas de seca, que ocorrem devido a 
fenomenos climaticos como el Nino, quando a variabilidade de alimento diminui sensivelmente. 
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2.2.2 Resultados 

As tabelas l2~Tl e 12.21 mostram os valores dos parametros do modelo utilizados nas simulag5es. 



Descrigao das grandezas para (macho e femea) 


Nome 


Valor 


Mutagao (na tira de bits) 


M 


1 (bit de 32 bits)* 


Dominancia 


D 


1 


Limite de doengas 


L 


1 


Numero de filhos 


NF 


10 


Made minima de reprodugao 


R 


9 



Tabela 2.1: Parametros relativos a parte estruturada por idade. 



Grandezas (para machos e femeas) 


Nome 


Valor 


Dominancia 


DF 


16* 


Mutacao do fenotipo num bit 


MF 


0.8 


Grau de competigao 


X 


0.5* 


Fenotipos intermediaries 


k 


16* 


Fenotipos extremos 


k 


< 16 e > 16* 



Tabela 2.2: Parametros relativos a parte nao estruturada por idade. Sao considerados fenotipos 
intermedidrios apenas aqueles com k = 16. *Grandezas que nao mudam neste capitulo. 

A distribuigao de fenotipos da populagao no primeiro periodo de tempo corresponde a uma 
gaussiana centrada em k = 16 (ver figura l2~T| quadrados). Este comportamento unimodal e 
estatico e deve-se ao numero de posigoes dominantes do fenotipo (DF = 16) e ao grau de 
competigao (X = 0.5). Neste intervalo de tempo nao existe uma selegao do fenotipo em fungao 
da fonte de alimento. 

A distribuigao de fenotipos gerada no segundo periodo de tempo pode apresentar-se de duas 
formas diferentes, para um mesmo conjunto de parametros (tabela l2~2~j) . dependendo apenas da 
semente aleatoria inicial: 

i) Ela pode permanecer em um estado estacionario oscilatorio entre duas distribuigdes bimodais, 
represent adas, respectivamente, pelos cfrculos e triangulos da figura 12.11 Esta variabilidade 
genetica e chamada de polimorfismo, termo que foi introduzido por E.B. Ford [2] em 1940 e esta 
associado a variabilidade genetica mediante selegao. Indica que a populagao de passaros como 
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Figura 2.1: Distribuigdo de fenotipos da populagdo. No primero periodo de tempo, antes da 
ecologia mudar, a distribuigdo e unimodal (quadrados); no segundo periodo de tempo (t > 40 x 
lO 3 ^ aparece um estado estaciondrio no qual a distribuigdo oscila entre aquela representada pelos 
circulos e aquela representada pelos tridngulos. As medias dos fenotipos extremos permanecem 
centradas em 6 e 26, respectivamente. 




Figura 2.2: Dindmica da populagdo de herbivoros no caso do estado estaciondrio oscilatorio da 
figura anterior. Lado esquerdo: pode-se observar o comportamento dos tamanhos das populagoes 
com fenotipos extremos, u pequenos" versus "grandes" ', durante toda a simulagao. O atrator 
semelhante a um borrao nesta figura corresponde aos tamanhos das populagoes dos extremos 
antes da ecologia mudar. Lado direito: pode-se observar como o tamanho de cada populagdo 
dos extremos varia com o tempo, depois da mudanga na ecologia. 
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urn todo carrega tanto os genes para tamanhos de bicos grandes quanto para pequenos. No caso 
da capacidade de alimento variar novamente, voltando a distribuigao inicial, o mesmo ocorrera 
com a distribuigao de bicos (voltara a ser unimodal), ja que nada impede o acasalamento 
entre bicos extremos oposto, gerando filhotes de bico medio. Observe que a existencia de um 
polimorfismo e uma condigao necessaria para que uma especiagao venha a ocorrer, caso alguma 
selegao sexual seja introduzida para provocar o isolamento reprodutivo entre fenotipos extremos 
opostos. 

Em ambas as distribuigoes bimodais da figura 12.11 as populagoes com fenotipos extremos 
sao maiores do que a populagao de fenotipo intermediario (k = 16). Contudo, ora a populagao 
de bicos grandes e maior que a de bicos pequenos, ora e o contrario. Este comportamento 
oscilatorio entre as populagSes dos extremos tern periodo igual a 10 passos de simulagao e 
pode ser visto na figura ()2.2}) . Observe que este periodo corresponde tambem ao numero de 
intervalos de tempo (passos) desde o nascimento (idade zero) ate a idade minima de reprodugao 
da populagao, R = 9. 

ii) O polimorfismo pode tambem desaparecer: a distribuigao de fenotipos evolui para um estado 
estacionario sem variabilidade, onde a distribuigao e tambem unimodal mas de largura menor 
que a inicial (ver figura l2~^j) . Neste caso a populagao dos extremos permanece constante e 
pequena, e a populagao do meio domina, apesar da escassez da fonte de alimento para estes 
fenotipos. 

As distribuig5es finais em (i) ou (ii), antes de serem atingidas, passam por um estado 
oscilatorio instdvel no qual a populagao de fenotipo intermediario chega a ter um tamanho 
comparavel aos das populagoes com fenotipos extremos (ver escala da figura l2~^| triangulos). O 
tempo que a distribuigao passa neste estado instavel pode, com uma pequena probabilidade, ser 
bastante longo, permitindo que se observe as mesmas oscilagoes nos tamanhos das populagoes 
dos extremos presentes no caso (i), mas cujas amplitudes variam no tempo, como mostra o lado 
direito da figura 12.51 

Nas simulag5es pode-se favorecer o aparecimento do polimorfismo oscilatorio estavel mudando- 
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Figura 2.3: Distribuigao de fendtipos no caso estaciondrio sem variabilidade. Lado esquerdo: 
distribuigoes antes e depois da ecologia mudar. Lado direito: tamanhos das populacoes com 
fendtipos extremos durante toda a simulagao. 



se o valor da taxa de mutagao fenotipica MF (veja tabela l2~3*j) . Para valores pequenos desta 
mutagao, por exemplo para MF = 0.01, de 10 simulagoes feitas com sementes aleatorias difer- 
entes obtivemos 10 estados estacionarios deste tipo. Nestes casos a populagao final e grande, 
ja que as medias dos fendtipos se concentram em valores proximos de e 32, para os quais a 
fonte de comida e maxima. 



MF 


Media dos fendtipos extremos 


Estabilidade 


1.0 


7 e 25 


7 oscilatorios e 3 unimodais 


0.8 


6 e 26 


2 oscilatorios, 7 unimodais e 1 instavel 


0.1 


1 e 31 


10 oscilatorios 


0.01 


e 32 


10 oscilatorios 



Tabela 2.3: Resultados de 10 simulagoes para diferentes taxas de mutagao do fenotipo. 



Uma outra forma de se favorecer o aparecimento de um polimorfismo estavel e aumentar 
o grau de competigao para o fenotipo do meio, por exemplo trocando X = 0.5 por X = 1.0. 
Neste caso os fendtipos intermediaries ficam de tal forma prejudicados que impedem que a 
distribuigao unimodal sem variabilidade (fig. I2.3J) seja atingida. 
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Figura 2.4: Polimorfismo instdvel no qual a distribuigdo de fenotipos oscila entre a curva preta 
(tridngulos) e a curva azul (circulos), para finalmente se tornar unimodal e estdvel como na 
figura\K3 




Figura 2.5: Lado direito: Tamanhos das populagdes de fenotipos extremos no estado em que 
o polimorfismo e meta-estdvel, medidos durante inumeros passos apos a ecologia mudar. Lado 
esquerdo: oscilagoes nos tamanhos das populagdes dos extremos, cujas amplitudes variam no 
tempo. 
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2.2.3 Selecao sexual e Especiacao 

Como mencionado no capitulo anterior, a selegao sexual e definida pela femea, pois e ela que 
escolhe o macho. As femeas com "gene seletivo" igual a 1 decidem, dentre NM machos, o 
melhor parceiro para acasalar; as femeas com "gene seletivo" igual a acasalam aleatoriamente. 
E importante lembrar que as femeas herdam este gene da mae com uma taxa de mutagao 
reversfvel, MS. 

Iniciamos as simulagoes com todas as femeas nao seletivas (bit de seletividade igual a zero 
para todas). Para qualquer valor de MF, MS = 0.001, NM = 6 e com grau de competigao X = 
0.5, obtem-se uma distribuigao estavel, com as duas populag5es de fenotipos extremos isoladas 
reprodutivamente: pode-se observar na figura 12.61 que a populagao de fenotipo intermediario 
desaparece completamente (lado esquerdo), ja que 95% da populagao de femeas se torna seletiva 
(lado direito). Esta estabilidade e obtida tanto com a regra com diregao como com a regra da 
diferenga para o acasalamento, descritas no final do capitulo anterior. Note que o numero 
de femeas seletivas aumenta rapidamente (entre 20 e 60 geragoes), como pode ser visto pelo 
detalhe da figura 12*1)1 lado direito. 
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Figura 2.6: Lado esquerdo: Distribuigao de fenotipos da populagao com selegao ecologica e 
sexual. Lado direito: densidade de femeas seletivas para duas regras diferentes de acasalamento. 
O detalhe ("inset") mostra a diferenga de velocidades na obtengdo do isolamento reprodutivo, 
dependendo da regra. 



Comparando-se as figuras 12.11 e 12.61 pode-se concluir que a populagao de fenotipo interme- 
diario da figura l2~T1 era formada apenas por bebes resultantes do acasalamento entre fenotipos 
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extremos diferentes. Voltaremos a este ponto na ultima segao. Aumentando-se o valor da taxa 
de mutagao MS, a densidade de seletivos diminui, mas os valores medios da distribuigao nao 
mudam e nem tampouco sua estabilidade. 

Ressaltamos ainda que tambem realizamos simulagSes onde os fenotipos intermediaries cor- 
respondiam a uma faixa de valores, e nao apenas a k — 16. Desde que estes fenotipos estivessem 
dentro do intervalo n% = 10 < k < ri2 = 22, os resultados obtidos foram qualitativamente os 
mesmos. Contudo, para ri\ < 10 e ri2 > 22 a probabilidade de se obter a especiagao diminui, 
tao mais quanto mais proximos os valores de n x e n<i estiverem dos extremos. 

2.3 Especiagao para predadores 

No processo de especiagao anterior a especie herbivora sofre o efeito direto da mudanga na 
ecologia, isto e, da mudanga na distribuigao de plantas. Na natureza existem tambem especies 
que sofrem o efeito indireto de uma selegao ecologica, como por exemplo, as especies predadoras 
que se alimentam do herbivoro que se alimenta destas plantas. Nesta segao nos estudaremos o 
processo de especiagao simpatrica neste tipo de cadeia alimentar, onde por exemplo, passaros 
(predadores) se alimentam de insetos (herbfvoros) , os quais se alimentam de plantas. 

Os insetos sao representados por seus codigos geneticos mas nao sofrem selegao sexual; 
alimentam-se das plantas cuja capacidade de sustentagao, F(k,t), e dada pela equagao ()2.1|) . 
A dinamica desta populagao ja foi analisada na segao anterior e as possfveis distribuigSes de 
fenotipos destes hervfboros, p h (k,t), sao as mesmas apresentadas nas figuras l2~Tl 12.31 e 12.41 A 
tabela l2~^l mostra os parametros do modelo utilizados nas simulag5es desta segao. Observe que 
os parametros relativos aos herbi'voros nesta tabela sao os mesmos utilizados na segao anterior 
(tabela Ej) . 

A capacidade de sustentagao da fonte de alimento dos predadores, que tambem sao repre- 
sentados por seus codigos geneticos mas que sofrem selegao sexual, e dada por: 

F p {k,t) = 2 xlO 6 x p h (k,t), (2.2) 

onde /ceo fenotipo dos herbivoros e t, o numero de passos da simulagao. Esta equagao liga di- 
retamente o comportamento da populagao de predadores ao da populagao de herbivoros. Desta 

32 



forma, tanto a competigao por alimento entre os herbfvoros como aquela entre os predadores e 
dada pela equagao ll.2[ mudando apenas a capacidade de sustentagao da fonte de alimento de 
cada especie. 

2.3.1 Resultados 



Grandeza 


Prcdador 


Herbfvoro 


M 


1 


1 


D 


10 


1 


L 


10 


1 


R 


10 


9 


NF 


2 


10 


DF 


16 


16 


MF 


0.5 


0.8 


MS 


0.001 


0.0 


X 


0.5 


0.5 



Tabela 2.4: Pardmetros do modelo para as duas especies, relativos tanto aos genomas 
cronologicos (parte superior da tabela) quanto aos genomas sera estrutura de idade (parte in- 
ferior da tabela). Os fendtipos intermedidrios correspondem aqueles com k = 16, tanto para 
herbivoros como para predadores. 



A distribuigao de fenotipos da populagao de predadores pode ser vista na figura l2~71 antes 
e depois da distribuigao de plantas mudar. Note que apesar da taxa de mutagao dos fenotipos 
dos predadores [MF = 0.5) ser diferente daquela dos herbfvoros (MF = 0.8), as distribuigoes 
dos fenotipos extremos dos predadores estao centradas em 6 e 26, isto e, nas mesmas posigoes 
que as populagoes extremas de herbfvoros (veja lado esquerdo da figura l2~T|) . 

A especiagao simpatrica dos predadores so ocorre quando o o polimorfismo oscilatorio dos 
herbfvoros e atingido, mas a recfproca nao e verdadeira, isto e, nem sempre o plimorfismo dos 
herbfvoros leva a especiagao dos predadores. Quando esta especiagao ocorre, a populagao de 
predadores com fenotipos extremos apresenta oscilagoes com perfodo igual ao das oscilagoes da 
populagao de herbfvoros tambem com fenotipos extremos 1 , independentemente de qual seja a 
idade minima de reprodugao dos predadores (ver tabela 

1 Lembre-se que o perfodo das oscilagoes da populagao de herbfvoros depende da idade minima de reprodugao 
dos mesmos - ver figura l2~2l 
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Figura 2.7: Distribuigao de fendtipos da populagao de predadores, que estd sendo afetada pela 
mudanga subita na fonte de alimentgdo dos herbworos. 
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Figura 2.8: La<io esquerdo: Populagao de herbworos com fendtipos extremos grandes em fungao 
da populagao de predadores tambem com fendtipos extremos grandes, durante toda uma sim- 
ulagao onde ocorreu a especiagdo dos predadores. Lado direito: oscilagdes nas populagdes de 
fendtipos extremos, tanto de herbworos quanto de predadores. 
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Figura 2.9: "Velocidades" de especiagdoem cadeias alimentares de duas (curva superior, ver- 
melha) e de tfes especies (curva inferior, azul). 

Como pode-se ver no lado direito da figura 12.81 existe uma defasagem de 3 passos entre 
as oscilagSes da populagao de herbfvoros com fenotipos extremos e aquelas da populagao de 
predadores com fenotipos extremos, defasagem esta que nao muda com nehum dos parametros 
da tabelal2~4l 

No caso em que a mutagao no fenotipo dos herbfvoros vale MF = 1.0 e a mutagao no 
fenotipo dos predadores vale MF = 0.01, obtem-se casos onde a populagao de herbfvoros 
apresenta um estado estacionario com polimorfismo oscilatorio, mas a populagao de predadores 
nao especia, isto e, apresenta so uma caraterfstica fenotipica, que tanto pode ser a de um dos 
extremos (k < 16) como a do outro (k > 16). 

A figura 12.91 mostra que numa cadeia alimentar de tres especies (plantas, herbivores e 
predadores) a especiagao e mais lenta do que numa cadeia alimentar de duas especies (plantas 
e herbfvoros), como era de se esperar. 

2.4 Perspectivas 

O polimorfismo oscilatorio instavel (figuras 12.41 e 12.5)1 . diffcil de ser observado, apresenta 
caraterfsticas de um estado meta-estavel (tempo de relaxagao grande). Alem disso, mantem 
os tres tipos de fenotipos (dois extremos e um intermediario) com populagoes de tamanhos 
comparaveis, pois neste caso a populagao intermediaria nao e composta apenas por bebes: 
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apresenta a mesma curva de mortalidade das populagoes com fenotipos extremos. A figura 
12.101 compara as distribuigoes por idade das populagoes intermedidrias, nos dois casos em que 
o polimorfismo aparece, o estavel e o instavel. Esta distribuigao tern a mesma forma daquela 
usualmente obtida com o modelo Penna e mostra que no caso estavel, de fato a quase totalidade 
dos hfbridos sao bebes resultantes do acasalamento entre fenotipos extremos, bebes estes que 
nao chgarao vivos ate a idade minima de reprodugao, R = 9. Ja no caso instavel, os hfbridos 
conseguem sobreviver o suficiente para reproduzir durante alguns perfodos. 

Seria interessante compreender que parametros interferem no tempo de vida deste estado 
meta-estavel e verificar se existe alguma faixa de parametros para a qual este polimorfismo se 
torna estavel. 
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Figura 2.10: Lado esquerdo: Distrubuigoes medias por idade das populagoes de fenotipo in- 
termedidrio, medidas nos ultimos 250 passos e normalizadas (divididas pelos seus respectivos 
valores mdximos). Pode-se observar que a longevidade dos hibridos e bem maior no caso instavel 
do que no estavel, em R = 9. Lado direito: as mesmas distribuigoes em escala semi-logaritmica 
mostram mais claramente a diferenga entre os dois casos. 



Os fato de que as oscilagSes na populagao de predadores (cadeia alimentar de tres especies) 
tenham perfodo igual ao das oscilag5es na populagao de herbfvoros, parece ser devido a relagao 
direta entre estas duas especies, dada pela equagao ()2.2j) . Contudo, nao conseguimos identificar 
o(s) parametro(s) responsavel(eis) pela diferenga de fase de 3 passos entre estas oscilagoes. 
Para entender melhor este comportamento gostarfamos de trabalhar numa cadeia alimentar de 
3 especies, na qual os predadores se alimentarao tanto dos herbfvoros quanto das plantas cuja 
distribuigao muda repentinamente, estando, desta forma, submetidos a uma selegao ecologica 
direta e a outra, indireta. 
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Capitulo 3 

Modelo com dois tragos fenotipicos 



Apresentaremos agora urn modelo no qual os individuos tern urn trago relacionado a selegao 
ecologica e urn outro trago, independente, relativo a selegao sexual. A inspiragao biologica 
surgiu a partir das farmlias de ciclideos encontradas na Nicaragua e na Africa e que vein 
atraindo grande interesse, em fungao dos dados estatfsticos coletados acerca destas fanrilias. 

Atraves de simulag5es primeiro estudamos a distribuigao dos fenotipos ligados a selegao sex- 
ual, numa ecologia onde a capacidade de sustentagao do ambiente independe do trago ecologico. 
objetivo e o de compreender melhor as consequencias da selegao sexual, independentemente 
da ecologia. A seguir, introduzimos a ecologia dependente do trago ecologico a fim de estu- 
dar o processo de especiagao simpatrica nestes peixes. Analisamos os resultados do processo 
para diferentes valores das taxas de mutagao, com o objetivo de comparar qualitativamente os 
resultados do modelo com aqueles encontrados experimentalmente nos ciclideos da Nicaragua. 

3.1 Motivacao Biologica 

Os ciclideos parecem especiar rapidamente, pois analises das 300 especies encontradas no lago 
Vitoria, no leste da Africa, indicam a existencia de um mesmo ancestral comum, que teria vivido 
ha apenas 12.400 anos atras j2Hj. Os estudos geneticos da especie Amphilophus citrinellum da 
Nicaragua jTT] revelaram que estes peixes apresentam distintos tipos de maxilares, que sao a 
chave da sua adaptagao nas diversas ecologias existentes. Um tipo de maxilar e adequado a 
uma dieta de caracois (molariform) e um outro tipo e adequado a ingestao de presas moles 
(papilliform). 
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Esperava-se que a especiagao simpatrica dos ciclideos envolvesse grande divergencia no trago 
ecologico e que houvesse uma forte correlagao entre este trago e a cor dos peixes, que e o trago 
ligado ao acasalamento. Contudo, nao e assim. As estatisticas mostram que existe uma forte 
diferenciagao na cor (entre os tipos "gold" - dourado e "normal" - listrado de preto e branco), 
mas que nao existe grande diferenciagao na morfologia dos maxilares, isto e, dentre os peixes 
de uma dada cor os dois tipos de maxilares sao encontrados, ao inves de apenas um tipo de 
maxilar para cada cor. De fato, a correlagao entre os dois tragos mostrou-se pequena (R = 0.48). 
Voltaremos a estas observagoes no final deste capftulo. 

Dois importantes modelos teoricos (do estilo "campo medio") foram desenvolvidos con- 
comitantemente [2E1 ETJ , para explicar a especiagao nos ciclideos. Ambos obtem a especiagao 
envolvendo dois tragos independentes, agregada a uma forte correlagao entre estes dois tragos 
[26 . Isto e, no caso dos ciclideos, os modelos prevem que a quase totalidade dos peixes doura- 
dos se alimentaria de caracois e os listrados, de presas moles (ou vice- versa). No trabalho |26- 
os autores admitem uma mudanga brusca tanto na ecologia como na competigao, para obter a 
especiagao simpatrica. A competigao so se torna intra-especifica quando a ecologia muda: antes 
disto, todos disputam igualmente com todos pela fonte de alimento. Esta mesma estrategia foi 
utilizada com sucesso nas simulagSes apresentadas em [20] e tambem nas nossas apresentadas 
em |2HJ [2H], na obtengao da especiagao. Contudo, no modelo adotado em [27] e considera- 
do mais realista |3U] , nao existem tais mudangas: tanto a ecologia quanto a competigao por 
alimento sao mantidas fixas e, ainda assim, a especiagao e obtida. 

Como visto nos capftulos anteriores, nos modelos que apresentamos nesta tese a competigao 
e a mesma desde o inicio da simulagao e dada pela equagao 11.21 com X = 0, 5. Contudo, 
a ecologia muda, de uma onde a abundancia de alimento e a mesma para qualquer fenotipo 
para outra que desfavorece fortemente os fenotipos intermediaries. Com o objetivo de suavizar 
tal desfavorecimento, adotaremos para o segundo intervalo de tempo a nova capacidade de 
sustentagao do ambiente, ilustrada na figura 13.11 e dada pela seguinte equagao: 



F(k,t) = 5 x 10 5 x 



1 

1.0 - 0.8 x exp(-(W - k) 2 /6A) 



< t < 4 x 10 3 , 

4 x 10 3 < t < 5 x 10 4 



(3.1) 
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Figura 3.1: Comparagao das capacidades de sustentagao do ambiente, no segundo intervalo 
de tempo (t > 4 x lO 3 ^, utilizadas no capitulo anterior (linha solida) e neste capitulo (linha 
pontilhada). 

Em nossas simulag5es os "ciclfdeos" tern a morfologia do maxilar representada por urn par 
de tiras de bits (fenotipo k) e a cor e tambem representada por urn par de tiras, independente 
do primeiro (fenotipo k'). 

3.2 Distribuigao do traco sexual numa ecologia fixa 

Como dito no initio deste capitulo, estudaremos primeiro o papel da selegao sexual sozinha, 
sem qualquer mudanga na ecologia. A seguir mostramos os resultados das simulagoes feitas 
com a mesma capacidade de sustentagao fixa anterior (t < 4 x 10 3 ), mas com individuos que 
tern agora mais uma caracteristica fenotipica, k', ligada a selegao sexual. Lembramos que as 
simulag5es comegam com toda a populagao de femeas nao seletiva. A tabela 13.11 mostra os 
parametros utilizados, onde MS e a taxa de mutagao no carater da seletividade. 

O comportamento da distribuigao dos fenotipos k nao muda em relagao aos resultados 
encontrados no capitulo anterior (ver figura I7!d1 lado esquerdo), isto e, a distribuigao novamente 
e uma gaussiana centrada em 16 e independente de MF. Ja a distribuigao dos fenotipos k! 
depende dos valores de MFS, MS e NM, como veremos a seguir. 
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Descrigao das grandezas para (macho e femea) 


Nome 


Valor 


Mutacao (na tira de bits) 


M 


1* 


Dominancia 


D 


3* 


Limitc dc doengas mortais 


L 


3* 


Numero de filhos 


NF 


2* 


Idade minima de reprodugao 


R 


10 


Dominancia - fenotipo k 


DF 


16* 


Dominancia - fenotipo k' 


DFS 


16* 


Mutacao por bit - fenotipo k 


MF 


0.1 


Mutacao por bit - fenotipo k! 


MFS 


0.5 


Fenotipos intermediaries 


kek' 


16* 


Mutacao na seletividade 


MS 


0.001 


Numero de machos dispomveis 


NM 


6 



Tabela 3.1: Valores dos pardmetros usados nas simulagdes. Os pardmetros relativos a parcela 
estruturada por idade do genoma estdo na parte superior. *Grandezas que nao raudara neste 
capitulo. 




Figura 3.2: Distribuigao do fenotipo k' ligado a selegao sexual. No inicio das simulagdes 
(t = 200 passos) a distribuigao fica centrada em 16 (circulos). O valor medio final se desvia de 
16 (quadrados), tanto para a direita quanto para a esquerda, dependendo da semente aleatoria 
inicial. Quanto menor o valor de MFS , mais a distribuigao se aproxima de um dos extremos 
(tridngulos) . 
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No imcio da simulagao muitas femeas ainda sao nao seletivas e por esta razao a distribuigao 
fica centrada em 16 - figura l3~2"l cfrculos. A medida que as femeas se tornam seletivas, uma das 
duas carateristicas fenotipicas extremas se fixa aleatoriamente - figuraE21 quadrados cheios ou 
vazios. Este desvio da distribuigao para um fenotipo extremo acontece tanto com a regra da 
diferenga para o acasalamento quanto com a regra com forte diregao. Observe que na regra com 
forte diregao a femea seletiva escolhe o macho mais ao extremo de sua propria cor, mas mesmo 
assim a distribuigao final nao apresenta uma divisao mas sim, um desvio. valor medio da 
distribuigao final e o tempo que esta leva para ser atingida depende dos valores das taxas de 
mutagao, MFS e MS, e do mimero de escolhas dispomveis na hora do acasalamento, NM. 
Em comparagao com os valores da tabela EIH se o valor de MS for muito pequeno, o valor 
medio fica em 16 por um tempo maior; se o mimero de escolhas NM for muito grande, o valor 
medio vai para um dos extremos (0 ou 32). 

3.3 Especiagao 

Como acabamos de ver, numa ecologia (ou capacidade de sustentagao) fixa e independente 
do fenotipo k e com um grau de competigao medio (X = 0,5), a selegao sexual nao e capaz 
de levar a especiagao. Introduziremos entao a mudanga na ecologia dada pela equagao 13.11 e 
observaremos como este cenario se modifica. 

Antes de apresentar os resultados precisamos definir o que chamaremos de especiagao, neste 
caso em que os tragos sao independentes. Observe-se que esta definigao e bastante problematica 
para os biologos Alguns acreditam que para haver especiagao e necessaria apenas a selegao 
sexual [3T] e neste caso, nao tern duvida quanto a especiagao dos ciclideos, pois estes apresentam 
uma clara divergencia na cor. Outros autores, entretanto, acreditam que a especiagao decorre 
da associagao entre uma selegao ecologica disruptiva e uma selegao sexual [2E1IH2]- Para estes 
autores, os ciclideos estao ainda em processo de especiagao, que so se concretizara quando 
os maxilares se dividirem, da mesma forma que a cor. Nos, contudo, temos uma vantagem 
com relagao ao mundo real: podemos medir exatamente a densidade de seletivos na populagao 
e saber quando uma parcela da mesma se torna reprodutivamente isolada, isto e, uma nova 
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especie. Assim, sempre que a densidade de seletivos final for proxima de 1 e constante no tempo 
(sem grandes flutuagoes), diremos que houve especiagao. 

Alem das distribuig5es dos fenotipos k e k', mostraremos tambem a correlagao, R(t), entre 
estes dois tragos. Ela e dada por: 

. > < kk! > - < k >< k! > 
R(t) = 

onde < k! >= Y^i=\^ K/P°p(t) e o valor medio da distribuigao de cores da populagao total no 
tempo t, Pop(t), e (Tfc/ e o desvio padrao deste valor medio. O mesmo tipo de calculo e feito 
para a morfologia k do maxilar e para o produto kk'. 

3.3.1 Resultados similares aos de Kondrashov |26j 

Na tabela IT21 mostramos os resultados obtidos em 10 simulagoes distintas, usando os valores 
das taxas de mutagao MF e MFS tambem mostrados na tabela. Em todos os casos a taxa 
de mutagao do gene da seletividade vale MS = 0.001. Os numeros entre parenteses correspon- 
dent a media, nas 10 simulag5es, dos valores medios finais da correlagao R e da densidade de 
femeas seletivas p s obtidos em cada simulagao. Os resultados aqui apresentados foram obti- 
dos utilizando-se a regra com forte diregao para o acasalamento, mas adotando-se a regra da 
diferenga os resultados sao os mesmos. 

Os casos que chamamos de especiagao correspondem aqueles em que as distribuigoes finais 
dos dois fenotipos se tornam bimodais estaveis e que, como dito anteriormente, apresentam uma 
densidade de femeas seletivas proxima de 1 e estdvel. Os casos IMaxilar/lCor correspondem 
aqueles em que ambas as distribuigoes se mantiveram unimodais depois da mudanga na ecologia. 
Nos casos assinalados como 2Maxilares/lCor apenas a distribuigao do fenotipo k se dividiu, 
tendo a distribuigao de k' se mantido unimodal. 

Foi escolhido um caso da tabela E21 onde ocorreu especiagao, para mostrar as ditribuigoes 
estacionarias dos tragos ecologico e sexual, a evolugao temporal das populagoes de cada fenotipo 
e os comportamentos da densidade de seletivos e da correlagao durante o processo de especiagao 
- figura EHfl (a) — > (e). Nos outros casos da tabela o comportamento das correspondentes 
grandezas sao os mesmos. No caso (f) da figura EPl a especiagao nao ocorreu e pode-se comparar 
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Figura 3.3: Especiagdo com MS = 0.001, MF = 0.1 e MPS = 0.5; a mudanga na ecologia 
se dd no passo t = 4 x 10 3 . (a) Distribuigdo do trago ecoldgico k. (b) Distribuigdo do trago 
sexual k'. Os circulos correspondem ao intervalo de tempo inicial (t < 4 x lO 3 ^ e os quadrados 
ao final da simulagao. (c) Evolugao temporal das populagoes com fendtipos k e (d) evolugao 
das populagoes com fendtipos k' , durante o processo de especiagdo. e) Variagao temporal do 
valor absoluto da correlagao R (linha solida) e da densidade de seletivos p s (linha pontilhada). 
(f) Variagao temporal da correlagao e da densidade de seletivos num caso onde nao houve 
especiagdo (2Maxilares/lCor) . 
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MF 


MFS 


Especiagao 
R = 0.985, p s = 0.997 


IMaxilar/lCor (R, p s ) 


2Maxilares/lCor (R, p s ) 


0.1 


1.0 


9 





1 (0.85, 0.6) 


0.1 


0.5 


2 


1 (0, 0.6) 


7 (0.85, 0.6) 


0.1 


0.1 








10 (0.8, 0.6) 
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0.5 





1 (0, 0.6) 


9 (0.85, 0.6) 


0.5 


0.1 








10 (0.5, 0.5) 
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0.1 








10 (0, 0.5) 



Tabela 3.2: Resultados, no final das simulagoes, para MS = 0.001. Em todos os casos de 
especiagao as distribuigdes finais do tragos k e k' sao bimodais e podem ser vistas na figura 
\3.tA (a) e (b). O estado IMaxilar/lCor indica que no final da simulagao foram encontradas 
distribuigdes unimodais tanto para k como para k! , centradas nos extremos. 2Maxilares/lCor e 
o caso em que a distribuigao de k' e unimodal e centrada num extremo e a distribuigao de k e 
bimodal. Os numeros entre parenteses sao a media dos valores medios finais de R e p s obtidos 
nas 10 simulagoes. 



as flutuagoes, tanto em R como em p s , com aquelas mostradas em (e), onde houve especiagao. 

Observando a tabela 13.21 nota-se que a probabilidade de especiagao diminui a medida em 
que MF aumenta. A especiagao e mais provavel quando MF e pequeno porque a distribuigao 
bimodal dos fenotipos k apresenta os dois modos mais separados do que para valores maiores 
de MF (ver tabela l2~3~|) . Uma maior separagao destes dois modos provoca um aumento nas 
populagoes dos extremos (que sao as seletivas) e uma maior eliminagao da populagao com 
fenotipo intermediario (que acasala aleatoriamente). 

Uma vez fixado um valor de MF, a probabilidade de especiagao diminui a medida em que 
MFS diminui. Valores pequenos de MFS difilcultam a especiagao porque fazem com que a 
distribuigao de cores (fenotipo k'), antes da mudanga na ecologia, se desloque do centro para 
um valor proximo a um dos extremos, como ja mencionado anteriormente por ocasiao da fig l3.2l 
e como tambem pode ser visto na figura l3~3l -fb). Observando a fig J3.3l -fd) ve-se que de fato uma 
das cores, no caso k' < 16, praticamente desaparece da populagao no momento da mudanga 
na ecologia, para reaparecer cerca de 30 mil passos depois e se igualar a populagao que havia 
permanecido. Pela fig l3.3l -fc). vemos que o mesmo nao ocorre com a distribuigao dos maxilares, 
isto e, assim que se da a mudanga na ecologia a populagao com um certo tipo de maxilar (no 
caso k < 16) diminui, mas rapidamente volta a crescer e estabelece-se um polimorfismo no que 
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diz respeito a morfologia dos maxilares. Contudo, o processo de especiagao so se concretiza 
quando a populagao com k' < 16 reaparece, quando entao a correlagao entre os dois tragos e a 
densidade de femeas seletivas vao para 1. Convem ressaltar que apresentamos tanto na tabela 
13.21 quanto na fig ]3.3l -fe) o valor absoluto da correlagao, mas que esta tanto pode ser positiva 
quanto negativa. Ela e positiva quando os extremos das populagoes de maxilar e cor coincidem, 
isto e, quando todos os individuos com k > 16 tern tambem k' > 16 e todos com k < 16 tern 
k' < 16. A correlagao e negativa quando os extremos sao opostos {k > 16 com k! < 16 e k < 16 
com k' > 16), mas em ambos os casos toda a populagao com um tipo de maxilar tern uma dada 
cor e toda a populagao com o outro tipo de maxilar tern uma mesma cor, diferente da primeira. 

Tanto a evolugao temporal dos fenotipos k e k' durante o processo de especiagao como o 
comportamento da correlagao que se estabelece entre eles sao bastante semelhantes aos obtidos 
em [26 , apesar dos modelos serem distintos. Contudo, nao parecem descrever corretamente 
o que ocorre com os cich'deos, ja que nestes observa-se uma clara divisao na cor mas nao na 
morfologia dos maxilares e ainda, as correlagoes entre os dois tragos sao pequenas se comparadas 
tanto com as obtidas em j2H| como com as nossas. 

Os processos de especiagao com MFS = 1.0 sao mais rapidos do que com MFS = 0.5, mas 
todos apresentam as mesmas carateristicas dadas pela figura l3~3~l (a) — > (e). 

3.3.2 Resultados mais adequados aos ciclfdeos 

Nesta subsegao adotaremos exatamente o mesmo procedimento da subsegao anterior, mas agora 
para uma taxa de mutagao do gene da seletividade MS = 0.0001, isto e, dez vezes menor. 

Comparando os resultados desta subsegao (ver tabela 13.31 e figura 13. 4j) com os da subsegao 
anterior, pode-se observar que a diminuigao de MS leva aos seguintes efeitos: 

1) O numero de casos em que ocorre especiagao aumenta (comparar tabelas). 

2) A distribuigao de cores permanece centrada no meio (k' = 16) ate o momento em que se da 
a mudanga na ecologia ffig l3.4l -(b^. 

Os ftens 1) e 2) corroboram a conclusao de que quanto mais centrada estiver a distribuigao 
de cores antes da ecologia mudar, maior a probabilidade da populagao especiar. 
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Figura 3.4: Especiagao para MS = 0.0001, MF = 0.1 e MFS = 0.5. (a) Distribuigao dos 
fenotipos k - trago ecoldgico. (b) Distribuigao dos fenotipos k' - trago sexual, (c) Comporta- 
mento da populagao com fendtipo k e (d) com fendtipo k! , durante o processo de especiagao. (e) 
Variagao temporal do valor absoluto da correlagdo R (linha solida) e da densidade de seletivos 
p s (linha pontilhada). (f) Variagao temporal da correlagdo e da densidade de seletivos num 
caso onde nao ocorreu especiagao(2Maxilares/lCor). 
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MF 


MFS 


Especiagao 
R = 0.974, p s = 1.000 


IMaxilar/lCor (R, p s ) 


2Maxilares/lCor (R, p s ) 
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1.0 
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2 (0.9, 0.55) 
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9 (0, 0.5) 


1.0 


0.1 
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Tabela 3.3: Resultados, no final das simulagdes, para MS = 0.0001. Nos casos de especiagao 
as distribuigoes finals dos tragos k eh! sao bimodais e podem ser vistas na figura \^[ (a) e (b). 



3) A correlagao entre os tragos ecologico e sexual diminui ligeiramente (ver tabelas l3~2l e 13. 3j) . 

4) A divisao na distribuigao dos maxilares ocorre ao mesmo tempo que a da cor e nao antes, 
como no caso anterior. 

Os itens 3) e 4) parecem indicar que valores menores de MF sao mais apropriados para 
descrever os ciclfdeos. Veremos nas conclusoes deste capitulo que e possfvel melhorar um pouco 
mais estes resultados. 

Nas tabelas 13.21 e 13.31 nao foi encontrado nenhum caso de especiagao para MF = 1.0 e 
MFS = 0.1, mas estes podem ser obtidos utilizando-se valores altos do numero NM de machos 
que a femea tern a disposigao para escolher na hora do acasalamento. Quanto maior for NM, 
mais proximo do momento em que a ecologia se divide se da o crescimento subito da densidade 
de seletivos e consequentemente, mais rapida e a especiagao, como mostra a figura 1331 

3.4 Perspectivas 

Como pode ser notado, neste capitulo nao foi dito se os estados IMaxilar/lcor e 2Maxilares/lcor 
sao estaveis, pois aumentando-se o tempo de simulagao para 160 x 10 5 passos as distribuigoes 
dos dois tragos nao mudam, mas a densidade de seletivos ffutua entre 0.1 e 1.0, em intervalos 
variados com e sem sicronia, dependendo da semente aleatoria inicial. Gostariamos de estudar 
estes estados, analisando as series temporais da densidade de seletivos e das populagSes com 
metodos como o TISEAN (Nonlinear Time Series Analysis) jHS]- O objetivo e conhecer o 
tipo de estabilidade dos casos de especiagao das tabelas e o tipo de instabilidade dos estados 
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Passos da simulagao * (10 ) 

Figura 3.5: Correlagao em fungao do tempo. O tempo que a densidade de seletivos leva para 
se aproximar de 1 depende do numero de escolhas NM. Quanto mais exigente e a femea, mais 
rdpida e a especiagao. 




Figura 3.6: Especiagao para MS = 0.0001, MF = 0.5 e MFS = 0.5. Lado esquerdo: dis- 
tribuigao dos fenotipos k - trago ecoldgico. Lado direito: distribuigao dos fendtipos k' - trago 
sexual. 
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IMaxilar/lCor e 2Maxilares/lCor. 

Gostariamos ainda de continuar procurando, na regiao de MS pequeno, uma combinagao 
de parametros que resulte numa correlagao mais baixa e numa divisao da cor anterior a dos 
maxilares, como ocorre nos cich'deos. Na figura ETo1 por exemplo, mudamos apenas o valor de 
MF = 0.1 para MF = 0.5; com isto, ja obtivemos uma populagao de fenotipos intermediaries 
k (maxilares intermediarios) que nao desaparece completamente, embora as cores ja estejam 
totalmente separadas. Contudo, a correlagao (nao mostrada) continua proxima de 1. 
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Capitulo 4 

Transigao de fase na especiagao 
simpatrica 

Neste capitulo utilizaremos urn modelo no qual o trago fenotipico relacionado a selegao ecologica 
e o mesmo utilizado para a selegao sexual e que comanda a competigao entre os indivfduos. A 
motivagao biologica deste modelo se relaciona com o aspecto da irreversibilidade do processo 
de especiagao em fungao do isolamento reprodutivo. 

Inicialmente mostraremos o comportamento da distribuigao dos fenotipos para diferentes 
graus de competigao X. Veremos que para valores altos de X, isto e, quando os fenotipos 
intermediaries estao submetidos a uma forte competigao, e possiVel obter especiagao sem qual- 
quer mudanga na ecologia. A seguir, estudando o valor medio da densidade de seletivos em 
fungao do grau de competigao, neste ambiente ecologico que independe do tempo e do fenotipo, 
observaremos um comportamento caracterfstico de uma transigao de fase onde a densidade de 
seletivos corresponde ao parametro de ordem e o grau de competigao, ao parametro de controle. 

4.1 Motivagao Biologica 

A larva da mosca da maga, Rhagoletis pomonella, e um dos exemplos de uma possiVel especiagao 
simpatrica em curso |31j. A Rhagoletis acasala preferencialmente perto da fruta que e sua 
planta hospedeira; diferengas nas preferencias por determinado hospedeiro podem resultar num 
completo isolamento pre-acasalamento entre individuos da mesma especie. A maga Crataegus 
spp. e o hospedeiro nativo da Rhagoletis pomonella na America do Norte, mas nos anos 80, 
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uma populagao nova surgiu na maga domestica, Malus pumila. Os cromossomas desta nova 
populagao foram mapeados e comparados aos da primeira especie. Os dados indicam uma 
inversao do polimorfismo, ou seja, ha uma possibilidade das duas populagoes estarem voltando a 
fundir-se numa unica [33] . Este comportamento nas larvas da mosca levam a uma questao muito 
interessante a respeito da reversibilidade ou irreversibilidade do processo de especiagao. Este 
tipo de questao se enquadra perfeitamente em fenomenos do tipo transigao de fase. 

4.2 Especiagao independente da ecologia 



Descrigao das grandezas (macho e femea) 


Variavel 


Valor 


Mutagao (por tira dc bits) 


M 


1* 


Dominancia 


D 


3* 


Limite de doengas 


L 


3* 


Idade minima de reprodugao 


R 


10* 


Numcro dc filhos 


NF 


2* 


Dominancia 


DF 


16 


Mutagao do fenotipo por bit 


MF 


0.01* 


Mutagao na seletividade 


MS 


0.001* 


Fenotipos extremos 


m 


13 


Fenotipos extremos 


n 2 


19 


Numero de machos dispom'veis 


NM 


50 



Tabela 4.1: Parametros do modelo, relativos tanto aos genomas cronologicos (parte superior 
da tabela) quanto aos genomas sem estrutura de idade (parte inferior da tabela). Tempo de 
simulagao de 4 x 10 4 passos. *Grandezas que nao mudam neste capitulo. 



Nesta segao serao mostradas as distribuigoes dos fenotipos k da populagao para diferentes 
graus de competigao X, num ambiente ecologico onde a capacidade de sustentagao nao muda 
no tempo nem depende do valor de k, isto e, F(k,t) = F = 2 x 10 5 . Observe que neste caso os 
ingredientes de nosso modelo sao os mesmos do modelo de campo medio de Dieckmann e Dobeli 
|2Zj, que como mencionado no capitulo anterior (segao motivagao biologica), e considerado mais 
realista que o modelo tambem do tipo campo medio de Kondrashov e Kondrashov Desde o 
inicio da simulagao os individuos estao sujeitos as competig5es por alimento e acasalamento, em 
fungao de suas caraterfsticas fenotfpicas, k. A simulagao comega com todos os individuos nao 
seletivos, mas a populagao de femeas pode se tornar seletiva e neste caso acasalar de acordo com 
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a regra com forte diregao. E importante lembrar que em nossas simulagoes dividimos a populagao 
em tres grupos, de acordo com os valores de k. grupo Pi consiste dos individuos com 
< k < n\; o grupo do meio, P m , corresponde aos fenotipos intermediaries onde ri\ < k < n 2 
e o terceiro grupo, P 2 , corresponde aos fenotipos com n 2 < k < 32. O tipo de competigao entre 
os individuos de cada grupo e dado pela equacao ll.2l Os valores utilizados para ni, n 2 e demais 
parametros podem ser vistos na tabela 14.11 Observe que agora os fenotipos intermediaries 
correspondem a uma faixa de valores, ao inves de apenas a k = 16. 

4.2.1 Resultados 
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Figura 4.1: Distribuigdo dos fenotipos nos primeiros 200 passos de simulagdo e distribuigdo 
estdvel no final da simulagdo, apos 4 x 10 4 passos. Lado esquerdo: Competigao favordvel 
para a populagao com fenotipo intermedidrio (X = 0.0). Lado direito: Competigao altamente 
desfav or duel para a populagao com fenotipo intermedidrio (X = 1.0). 



Quando o parametro de competigao X = 0.0, os individuos com fenotipos intermediaries 
sao favorecidos pois nao competem com os grupos fenotipicos extremos ( equacao 11.2)1 . Neste 
caso a distribuigao de fenotipos corresponde a uma gaussiana estavel e independente da semente 
aleatoria inicial e que pode ser vista na figura 14.1) lado esquerdo. O comportamento da densi- 
dade de seletivos em fungao do tempo mostra que a populagao atinge um estado estacionario 
onde nao existem femeas que reproduzam seletivamente, ou seja, a populagao nao especia - ver 
figura l4~2) lado direito, curva inferior (azul). 

Contrariamente, quando introduzimos um forte grau de competigao X = 1.0 a populagao 
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com fenotipo intermediario e severamente desfavorecida, pois compete com as duas populagoes 
inteiras de fenotipos extremos. Novamente os resultados independem da semente aleatoria 
inicial e a distribuigao final dos fenotipos da populagao e estavel. Esta distribuigao estavel 
mostra que a maioria da populagao tern fenotipo k = ou k = 32, ou seja, os fenotipos 
intermediaries desaparecem como pode ser visto na figura 14.11 lado direito. O comportamento 
da densidade de seletivos em fungao do tempo mostra que a populagao atinge novamente um 
estado estacionario, mas que neste caso consiste apenas de femeas seletivas - ver figura l4~2"l lado 
direito, curva superior (vermelha). Vemos portanto que com X = 1.0 a populagao especia sem 
necessidade de uma mudanga ecologica, devido ao forte grau de competigao. 
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Figura 4.2: Lado esquerdo: Distribuigao dos fenotipos da populagao em tres diferentes in- 
stantes: apos 200 passos (circulos vermelhos), apos 3900 passos (triangulos pretos) e no final 
da simulagao, apos 4 x 10 4 passos (quadrados azuis). Nos tres casos X = 0.5. Lado dire- 
ito: Comparagao da densidade de seletivos em fungao do numero de passos de simulagao, para 
diferentes graus de competigao X. 



As distribuig5es dos fenotipos para um grau de competigao intermediario X = 0.5 podem 
ser vistas na figura l4*72| lado esquerdo, em tres diferentes instantes. A correspondente densidade 
de seletivos em fungao do tempo apresenta grandes flutuagoes, como mostra a figura l4~2"| lado 
direito, curva do meio (preta). 
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4.3 Transigao de fase 

Como acabamos de ver, para X = 0.0 a populagao nao especia, isto e, a densidade de seletivos 
e estavel e proxima de zero. Para urn grau de competigao alto, X = 1.0, a populagao especia 
com uma densidade de seletivos estavel e proxima de 1. Contudo, para X = 0.5 a densidade 
de seletivos flutua muito, o que sugere a existencia de uma transigao de fase entre um estado 
onde toda a populagao permanece nao seletiva e outro onde todas as femeas se tornam seletivas 
e ocorre a especiagao. Nesta segao estudaremos mais detalhadamente como se comporta a 
densidade de seletivos em fungao X, a fim de verificar se de fato e possivel afirmar a existencia 
e caracterizar esta transigao de fase. 
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Figura 4.3: Lado esquerdo: Comporiamento da densidade media de seletivos, obtida com 10 
simulagoes distintas, em fungao do grau de competigao X. Lado direito: Comportamento do 
desvio padrdo da densidade media de seletivos. 



Para cada valor de competigao X foram feitas 10 simulagoes mudando-se apenas a semente 
aleatoria; em cada simulagao foi calculado o valor medio da densidade de seletivos, p s , nos 
ultimos 10 4 passos. Calculamos entao a media destes 10 valores medios, cujo comportamento 
em fungao do grau de competigao pode ser visto na figura fOl lado esquerdo. No lado direito 
da figura pode-se observar o comportamento do desvio padrao desta media em fungao de X. 

Em todas as simulagoes anteriores a carateristica fenotipica de cada indivfduo era dada 
por um par de tiras de 32 bits cada. Reduzindo o tamanho das tiras de 32 para 16 bits, a 
dominancia fenotipica de DF = 16 para DF = 8 e a definigao dos fenotipos k intermediaries 
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de ni = 13 e n 2 = 19 para rii = 7 e n 2 = 9, observamos que os comportamentos de p s e o~(p s ) 
variam ligeiramente no lado X < 0.5, como mostra a figura l4~H 
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Figura 4.4: Comparando, para diferentes tamanhos de tiras de bits, os comportamentos da 
densidade de seletivos media (lado esquerdo) e de seu desvio padrao (lado direito), em fungao 
do grau de competigao X . 

Associando a densidade de seletivos ao parametro de ordem e o grau de competigao X ao 
parametro de controle, vemos que de fato a figura I4.3I indica a existencia de uma transigao de 
fase em torno de X c = 0.5, entre um estado onde e livre a troca de genes entre os individuos 
(p s = 0) e um outro estado no qual coexistem duas populagSes que nao acasalam entre si 
(p s = 1). Se esta e de fato uma transigao de fase, ela se da fora do equilfbrio (SEj, uma vez que 
os individuos da populagao do modelo de Penna original vivem num estado meta-estavel. 

A densidade de seletivos se comporta como o parametro de ordem pois muda de um valor 
zero para um diferente de zero a partir de X = X c e sua derivada, o desvio padrao de seu valor 
medio, mostra um comportamento singular em X c = 0.5. Esta singularidade no desvio padrao 
aumenta com o tamanho das tiras de bits usadas para caracterizar o fenotipo (ver figura l4~4l 
lado direito), indicando que o comprimento das tiras de bits faz o papel do tamanho do sistema 
numa transigao de fase usual. 

O grau de competigao critico X = 0.5 pode ser determinado analfticamente, como veremos 
a seguir. Observamos pelas simulagoes que as populagoes de fenotipos extremos tern o mesmo 
tamanho, Px(k,t) ~ P 2 (k,t), como era de se esperar, devido a simetria com que foi definida 
a competigao intra-especifica dada pela equagao 11.21 Assim, colocando Pi = P 2 = P e nesta 
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equagao e chamando V\ = V% — V e , obtemos que: 

V e (k,t) = (P e + P m )/F, 
V m (k,t) = (P m + X*2P e )/F. (4.1) 

Pela equagao acima vemos que se X = 0, V m < V e e portanto a probabilidade de sobre- 
vivencia dos fenotipos intermediaries e maior que a dos fenotipos extremos. J a para X — 1, 
ocorre o oposto. 

O valor de X para o qual a probabilidade de sobrevivencia ou, equivalentemente, a com- 
petigao entre todos os individuos e a mesma, pode ser imediatamente obtido impondo-se que 
V m (k, t) = V e (k,t) na equacao l4.1l e vale X c = 0.5. 




0.2 0.4 0.6 0.8 1 0.1 1 

X Log X 

Figura 4.5: Lado esquerdo: Densidade media de seletivos em fungdo de X para diferentes 
numeros de escolhas no acasalamento, NM. Lado direito: Desvio padrao da densidade media 
de seletivos em fungdo de X para diferentes valores de NM. 

Diminuindo o numero de machos que a femea tern para escolher ao acasalar de NM = 50 
para NM = 3, verificamos que o comportamento da densidade media de seletivos e de seu 
desvio padrao em fungao de X variam, conforme mostra a figura I431 

Observe que para valores pequenos de NM a transigao de fase e destruida. Desta forma, 
podemos associar o papel do numero de machos dispomveis para o acasalamento ao inverso do 
campo numa transigao de fase usual. 

Com o objetivo de caracterizar esta aparente transigao de fase medimos, para 100 simulagoes 
distintas no intervalo X — X c = X — 0.5 G [0.001,0.01], a densidade media de seletivos e seu 
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Figura 4.6: A inclinagao das retas foram obtidas a partir de 100 simulagoes distintas, ao inves 
de 10. Lado esquerdo: Comportamento, perto do ponto critico, da densidade media de seletivos. 
Lado direito: Comportamento, perto do ponto critico, do desvio padrao da densidade media de 
seletivos. 

desvio padrao. As leis de potencia encontradas podem ser vistas na figura H~6l A primeira, do 
lado esquerdo da figura, relaciona-se ao comportamento do parametro de ordem perto do ponto 
critico, ou seja, p s (X ph X c ) ~ \X — X C \P. A segunda, do lado direito, esta relacionada com o 
desvio padrao tambem proximo ao ponto critico, isto e, c Ps (X m X c ) ~ \X — X c | -7 . 

Sabemos que os resultados aqui apresentados nao sao conclusivos no que diz respeito a 
ordem da transigao. Para concluir se esta transigao e de primeira ou de segunda ordem e 
necessario, por exemplo, repetir as simulag5es que deram origem a figura l4~4l para tamanhos de 
tiras de bits maiores e verificar o que ocorre com o pico da suscetibilidade (desvio padrao da 
densidade media de seletivos). Se este pico divergir, a transigao e de segunda ordem, se nao, 
e de primeira. Esta perspectiva futura de trabalho nos parece de fundamental importancia, 
pois sinalizara se o processo de especiagao pode ou nao ser reversivel, apesar do isolamento 
reprodutivo inerente ao mesmo. Contudo, tal perspectiva demandara um tempo razoavelmente 
grande, visto que exigira modificagSes bastante complexas nos programas. 

Note que nos dois capftulos anteriores tmhamos ri\ = n 2 = 16 e neste capitulo tomamos 
n\ = 13 e n2 = 19. A razao de utilizarmos este intervalo foi porque percebemos que com 
estes valores os pontos das curvas presentes na figura l4~T)l apresentavam uma menor dispersao. 
Observamos ainda que se neste modelo forem usados valores altos de MF (MF pa 0.1 ou maior) 
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a especiagao obtida para X > 0.5 e perdida. 

Numa transigao de fase no equilfbrio jSZj, proximo ao ponto critico, acontecem mudangas 
drasticas no comportamento da maioria das grandezas caracteristicas do sistema, como na 
entropia, na entalpia, etc. Aqui tambem foram encontradas mudangas no comportamento de 
outras grandezas alem da densidade de seletivos, como por exemplo, na populagao total, na 
populagao com fenotipo intermediario, nas populagoes com fenotipos extremos, etc. Mas existe 
uma grandeza que se mantem aproximadamente constante durante todas as simulagoes que 
e P e + P m , onde P m e a populagao com fenotipos intermediaries e P e e qualquer uma das 
populagoes com fenotipos extremos, como mostra a figura I4~71 Esta constante foi conflrmada 
pela versao preliminar do modelo analitico que vem sendo desenvolvido ha bem pouco tempo 
principalmente pelo estudante de doutorado Veit Schwammle. 
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Figura 4.7: Numa transigao de fase muitas grandezas sofrem mudangas de comportamento 
proximas ao ponto critico, como ocorre por exemplo com a populagao de fenotipos intermediario s 
(circulos vermelhos) e com qualquer uma das populagoes de fenotipos extremos (quadrados pre- 
tos). Contudo, a soma destas duas populagoes, P m +P\ ou P m +P2, permanece aproximadamente 
constante (tridngulos azuis). 

Nesta versao preliminar utilizamos as seguintes equag5es diferenciais, do tipo presa-predador 
de Lotka-Volterra: 



A = aP 1 + bP m -c(P 1 + P m )P 1 



(4.2) 



p 



aP m + bP x - c(2XP 1 + P m )P„ 
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O primero termo destas equagoes representa o crescimento das populagoes em fungao dos 
nascimentos. Se so houvesse este termo, o crescimento seria expoencial (a > 0). O segundo 
termo representa o resultado das mutag5es e cruzamentos que fazem com que individuos de 
uma populagao gerem filhos com fenotipos correspondentes ao da outra populagao. Os ultimos 
termos destas equag5es correspondem aos fatores de Verhulst (equagao 11.2)1 das populagSes de 
acordo com seus respectivos fenotipos. Este modelo, apesar de muito simples, vem confirmando 
varios resultados obtidos com as simulagoes e sera um dos temas da tese do Veit. 

4.4 Perspectivas 

Conforme ja mencionado, foi demonstrada qualitativamente a existencia de uma transigao de 
fase no processo de especiagao simpatrica, mas nao foi encontrada a ordem desta transigao. 
Pretendemos nos aprofundar nesta questao, nao com o intuito de comparar nossos resultados 
com dados reais, mas com o objetivo de estudar a reversibilidade ou nao deste processo. Lem- 
bramos que o mimero de machos disponiveis, ligado ao isolamento reprodutivo, faz o papel do 
campo na transigao de fase encontrada. 

No que diz respeito a especiagao da mosca da maga, tal estudo nao tern apenas uma im- 
portancia teorica, mas tambem economica. Por exemplo, a fruticultura brasileira tern prejuizos 
da ordem de R$ 360 milhoes anuais (US$ 120 milhoes) por causa da mosca-da-fruta, que so 
se reproduz uma vez. O presidente da Biofabrica Moscamed Brasil disse que [38] a partir de 
setembro do proximo ano o Brasil produzira moscas estereis, que serao utilizadas no combate 
a mosca-da-fruta. 

Pretendemos tambem utilizar este modelo onde a capacidade de sustentagao e fixa e inde- 
pendente de k, para estudar as distribuigoes dos fenotipos para diferentes graus de competigao, 
mas agora considerando dois tragos fenotipicos independentes. 
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Capitulo 5 
Conclusoes 



Comparar os modelos que descrevem o processo de especiagao simpatrica no geral e ambiguo, 
pois como se viu nos capftulos 2 e 3, o modelo depende da especie que esta em processo de 
especiagao e quais sao as condigoes ecologicas pelas quais passa a dita especie. Por exemplo, 
dizer que o modelo com dois tragos fenotfpicos e mais realista que o modelo com apenas um 
trago e verdade para o caso dos cichdeos, pois eles apresentam duas carateristicas fenoripicas 
claramente distintas. Contudo, isto ja nao e verdade para os tentilhoes, pois a carateristica 
fenoti'pica ligada a competigao e a mesma ligada a selegao sexual, uma vez que o som, atraves 
do qual a femea escolhe o parceiro, esta diretamente relacionado a morfologia do bico. 

Uma das perguntas mais ti'picas em ciencia e se o modelo utilizado descreveu ou nao alguma 
carateristica especial real do fenomeno que foi estudado. No caso da especiagao dos tentilh5es, a 
bibliografia descreve com riqueza diversas carateristicas qualitativas destes passaros, mas ainda 
nao existem dados quantitativos suficientes. A seguir ressaltaremos os resultados qualitativos 
das simulagoes que coincidem com aqueles observados nos tentilhoes. O resultado qualitati- 
ve que mais chama a atengao e a rapida adaptagao da morfologia do bico a uma mudanga 
climatica abrupta (ver figura 12.61 "inset" ) e a existencia de um comportamento oscilatorio no 
aparecimento das populagoes com bicos extremos (ver figuras I2~T1 e 12 .2} lado direito). Note 
que na figura |U da introdugao, as oscilagSes no tamanho dos bicos dos tentilh5es estao dire- 
tamente relacionadas as oscilagSes climaticas. Por outro lado, nossos resultados sugerem que 
so uma mudanga abrupta da ecologia produz oscilag5es com periodos de uma geragao. Por 
esta razao pretendemos continuar procurando na bibliografia dos tentilhoes dados sobre estas 
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oscilagoes, com o objetivo de introduzir no modelo oscilagoes periodicas da ecologia e verificar 
se as oscilagoes das populagoes extremas seguem as da ecologia, ou se o periodo de uma geragao 
encontrado aqui se mantem. 

No caso dos ciclideos existem dados geneticos e uma descrigao detalhada da forma de acasala- 
mento e alimentagao destes peixes, como foi visto na introdugao e no capitulo 3. A principio 
foi dificil entender o processo de especiagao dos ciclideos, primeiro porque na bibliografia o 
numero de palavras tecnicas e muito grande e segundo porque as simulagoes com o modelo 
de dois tragos demandam muito trabalho, no sentido de que existem varios parametros para 
"sintonizar" a nm de se obter a especiagao. Como foi visto nas tabelas, 13.21 e I3.3[ a obtengao 
da especiagao nao so depende dos parametros mas tambem da semente aleatoria com que e 
iniciada a simulagao. Finalmente este tema ficou claro, pois agora percebemos que para repro- 
duzir a especiagao dos ciclideos e preciso obter primeiro uma separagao na morfologia da cor 
e posteriormente, uma separagao na morfologia dos maxilares. Nos obtivemos dois resultados 
diferentes: num deles a separagao dos maxilares acontece primeiro, da mesma forma que em 
|2*o] . e no outro os maxilares e as cores se separaram ao mesmo tempo. Pretendemos contin- 
uar procurando os efeitos "microscopicos" que produzam a separagao nas cores antes da dos 
maxilares. Para isto pretendemos estudar a estabilidade das distribuigSes dos fenotipos, cor 
e maxilar, como foi exposto nas perspectivas do capitulo 3 e suavizar a mudanga brusca na 
ecologia, ja que as mudangas climaticas detectadas nos lagos onde foram estudados estes peixes 
nao parecem ser muito relevantes. 

O resultado encontrado no capitulo 4, acerca do comportamento da densidade de seletivos 
em fungao da competigao, e bastante interessante do ponto de vista da compreensao dos modelos 
para especiagao simpatrica em geral. Por exemplo, existe uma grande diferenga entre os modelos 
(2Iij e j2Z|, ja que o modelo [2E1 precisa de uma mudanga ecologica radical enquanto o modelo j2Ij 
nao utiliza mudanga alguma e os dois obtem especiagao simpatrica. Agora que conhecemos o 
comportamento tipo transigao de fase na especiagao, compreendemos que o modelo [26] precisa 
da mudanga na ecologia devido a fraca competigao a qual os fenotipos intermediaries estao 
sujeitos, o que equivale a regiao X < 0.5 na figura l4~!?| lado esquerdo. Ja o modelo |2Zj, que 
utiliza uma forte competigao, equivale a regiao X > 0.5 da mesma figura, nao precisando de 
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qualquer mudanga radical na ecologia. 

comportamento do tipo transigao de fase na especiagao simpatrica tambem nos fez perce- 
ber o tipo de estudo que pode ser realizado na mosca da maga, como foi explicado no capftulo 
4. Alem do ponto de vista economico, o tenia da irreversibilidade na especiagao simpatrica 
e tambem interessante do ponto de vista teorico, pois a identificagao da irreversibilidade ou 
reversibilidade do processo da a oportunidade de se saber se o sistema esta dissipando energia 
e que grandeza esta fazendo o papel da energia neste modelo. 
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